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Radial Distribution Study of Vitreous Barium Borate 1 

Arthur Bienenstock, 1 ' Aaron S. Posner, and Stanley Block 

(January 6, L960) 

X-ray diffraction, radial distribution studies of a 20 percent barium oxide, SO percent 
B 2 < >:< glass have been performed using both the atomic and electronic distribution functions. 
From these distributions, the average barium-barium distance has been determined as 6.70 A. 
This distance is in good agreement with that predicted by Levin and Block on the basis of 
a structural interpretation of immiscibilit y data. 



1. Introduction 

Levin and Block [1] have attempted to explain 
immiscibility of alkali and alkali-earth metal oxides 

in borate and silicate glasses by means of crystal 
chemical principles. They considered two short- 
range structures proposed by Warren and Pincus [2] 
for the system at the limit of miscibility. For the 
particular system investigated here, BaO-B 2 3 , the 
type A struct lire, shown in figure 1, has the boron 
coordinated either triangularly or (etraliedrnlly by 
oxygen. Two bariums are bonded to the same oxygen 
with approximately a 180° bond angle. Using ionic 
radii, they calculate a separation of barium atoms 
of 5.50 A." 

The type B structure is shown in figure 2. Here, 
each boron is tetrahedrally coordinated by oxygen. 
The boron tetrahedron shares edges with ;i distorted 
barium cube. The predicted barium-barium sepa- 
ration in this case is 6.67 A. 
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Figure 1. Type A coordination. 

Boron can be either in tetrahedral or triangular coordination. The Ba-Ba 
distance is the sum of the Ba-0 distances. 
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With a simple model of the barium atoms forming 
a cubic lattice and a knowledge of the density and 
composition at the limit of miscibility, Levin and 
Block found that the average separation of barium 
atoms is 0.82 A. Hence, they predicted that the 
structure is type B. 

This paper describes the determination of the 
Ba-Ba distance by means of X-ray diffraction radial 
distribution techniques. It will be shown that the 
structure cannot be type A and agrees with type B. 

2. Experimental Procedures 

Measurements were made with a North American 
Philips X-ray dill'ractomoter, modified so that the 
volume defined by the cylindrical scat ten- shield 
could be evacuated. Intensities were determined at 
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Figure 2. Type B coordination. 

Boron must be in tetrahedral configuration. The boron and barium polyhcdra 
share edges. 
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each angular position by the average of two con- 
secutive measurements of the time taken for 6,400 
counts. The maximum difference between the two 
measurements was 3.5 percent. Readings were 
taken at intervals of 0.05 in s=4t sin 0/A in the 
range 0.3 to 3.5 and intervals of 0.1 in the range 
3.6 to 7.9 using CuKa radiation with an argon filled 
Geiger counter, and intervals of 0.1 in the range 6.1 
to 10.0 using MoK« radiation and a krypton filled 
Geiger counter. 

A warmup time of 24 lir was allowed, after which 
the unit was kept running until all the data for each 
radiation were obtained. The intensity at a standard 
position was checked periodically to guard against 
fluctuations in the incident intensity. 

Monochromatization of the reflected X-ray beam 
was accomplished by means of balanced filters [3]. 

The sample, supplied by G. W. Cleek and E. H. 
Hamilton of NBS, contained 20 mole percent BaO 
and 80 mole percent B 2 3 . The observed intensi- 
ties are shown in figure 3. 
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Figure 3. Diffraction pattern of vitreous barium borate. 

3. Theory 

Neglecting polarization and atomic constants, the 
amplitude of X-rays scattered from an electron dis- 
tribution p(R) is given by 



F(k)={ P (R) exv{-2irik-R)dR 



(1) 



where k is the reciprocal space vector. The inten- 
sity is given by 

I(k)=F(k)F*(k), 

where F*(k) is the complex conjugate of F(Jc). Then, 
using the convolution theorem [4], the Fourier trans- 
form of iQc) is given by 

P(£)= fp(r)p(r+^dr= (V(l)exp {2irik-R)dK (2) 



where, as in most X-ray diffraction transform tech- 
niques, we have assumed that the sample is of infinite 
extent. 

If the sample is such that I{k) is spherically sym- 
metric in fc-space, e.g., a powder or an amorphous 
material, the integral over the angular portions, ex- 
pressing k in a spherical polar coordinate system 

with the k z axis parallel to it, can be performed 
directly. In this case, 

P(R)= I{k)exp(2irikRcosa)k 2 smad(j>dadk 

Jo Jo Jo 

=2ir I (k) exp (2wi JcR cosa) k 2 sina da dk 

Jo Jo 

= (2//?) [ I (k) sin(2dfc#) k dk 
P (7?) = (2/ 4/iV) r s I(s) sin(sR) ds (3) 

where s = 2irk = 4Tr sin 0A. 

Here, P(R) is the three-dimensional Patterson 
function. Since, for a powder or an amorphous ma- 
terial, it is spherically symmetric, it may be written 
as P(R)/4:tR 2 , where P(7?) is the relative probability 
of finding two electrons separated by a distance R, 
and is called the electron radial distribution function. 
Then 

P (R)=(2R/tt) r s I(s) s'm(sR) ds. ( L) 

With the usual [5] considerations about the inabil- 
ity to make measurements of / for small s, and the 
subtraction of the contribution of each atom convo- 
luted with itself to P(R), this can be written 

P(R)-Po(R)=(2R/tc) f s{I(s)-2 m fJ}sm(sR)ds 

J Smin 

(5) 

where IZmjm is the weighted sum of the squares of the 
atomic scattering factors of the atoms in the sample 
and Pq(R) is the distribution obtained from a homo- 
geneous array of electrons with the same average 
electron density as the sample. 

If the usual assumption is made that the different 
atomic scattering factors can be expressed as mul- 
tiples of a normalized average scattering factor, / e , 
where normalization implies f e (s = 0) = \, then divi- 
sion of the integrand by the square of the normalized 
average scattering factor reduces the distribution to 
that of point atoms and yields the atomic radial 
distribution function commonly used [6]. This can 
be written 

J2 m X m ^R 2 (Pn-Po) = (2R/7r) f " H(8) sm(Rs)ds (6) 
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where i(s) = (/(«) — S^)//?, 2 m A„, is the effective 
number of electrons of atom m summed over atoms m, 
and Airli-p^ll! is the Dumber of atoms, each multi- 
plied by its effective number of electrons, between 
R and R+dRj from atom m. 

4. Calculations 

The coherent scattered intensity is obtained from 
the experimentally measured intensity by correction 
for polarization and incoherent scattering. The in- 
coherent scattering factor for barium was calculated 
using the equation [7] 



--Z-f/Z. 



(7) 



The coherent scattering factors for barium, oxygen, 
and boron as well as the incoherent scattering factors 
for boron and oxygen were obtained from standard 
references [8]. The incoherent scattering factors 
were corrected with the Breit-Dirac factor [7]. 

The density of the sample was 2.83 g/cm 3 and from 
this value and tin 4 molecular weight, the average 
electron density was calculated to be 0.79 elec- 
trons/A 3 . 

The data were placed on an absolute scale by 
setting 

S«J(«)*iA*=S { Z>il+^calc In. } &*• (8) 

After corrections for polarization and incoherent 
scattering had been applied, the functions s{I(s) 
— X m fni] and si(s) were computed from the experi- 
mental data. The integrals were replaced by sums 
of the form 

10.0 

S s v(.s') ain($R)A$ 

*=0.3 

which were performed on a high speed digital com- 
puter for 0.1 A spacings in It from R=0 to #=9.9. 
Multiplication by 2Ii/ir gave the right-hand sides of 
(5) and (6), and from these the radial distribution 
functions were calculated. 

5. Conclusions 

Figure 4 shows the electron radial distribution 
P(R)-P (R) obtained. The most significant features 
for the purpose of this study are the absence of a 
peak at 5.5 A and the large peak at 6.8 A. The 
peak at 6.8 A corresponds to the type B barium- 
barium separation prediction by Levin and Block. 
5.5 A is the barium-barium separation for the type A 
structure. 

Figure 5 shows the electron radial distribution 
P(R). The poor resolution doesnot allow for a cal- 
culation of the area under each peak. However, 
the areas may be calculated by determining the area 
between minimums of the peaks of the P(R)-P (R) 
distribution and adding to this area the integral of 
Po(R) over the range between the minimums. This 
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Figure 4. Electron radial distribution function, 4irR 2 [p(R)-p 
(R)]=P(R) — P (R), calculated from the data of figure 3. 
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R, Angstroms 



Figure 5. Electron radial distribution function, 47rR 2 p(R) = 
P(R), calculated from the data of figure 3. 



method suffers considerably in accuracy due to over- 
lap, but offers some means of calculating areas. For 
this reason, the areas of the first three peaks only 
were calculated. The results were 322 (electrons) 2 
for the 1.2 A peak, 1750 (electrons) 2 for the 2.55 A 
peak, and 2779 (electrons) 2 for the 3.1 A peak. 

A more complete picture of the type B structure 
is shown in figure 6. On the basis of this model, the 
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FrGURE 6. Schematic representation showing the coordination 
of oxygen atoms about the barium. 

The boron and barium configurations share edges. 



peaks in the radial distribution can be explained 
readily. The 1.2 A peak corresponds to the boron- 
oxygen distance. The predicted area of this peak is 

4X1.6X10X2X2=256(electrons) 2 

Here 4 is tin 4 expected number of oxygens surround- 
ing each boron, 1.6 is the number of borons in the 
basic unit used for scaling the intensities, 10 and 2 
are the number of electrons per ion associated with 
the oxygen and boron, respectively (these values are 
used throughout the paper), and the last factor of 2 
takes into account the fact that each peak appears 
twice in the distribution function. The observed 
area of the peak is more consistent with an oxygen 
coordination of 5 around boron. However, this co- 
ordination would be quite unusual and is inconsistent 
with the short B-0 distance. Hence it seems likely 
that some error in the distribution function obtained 
from the data has changed both the area and posi- 
tion of the peak. In the region of small R this func- 
tion is most susceptible to such errors and these 
errors will be most significant because of the small 
area of the peak. 

The 2.55 A peak corresponds to the oxygen-oxygen 
separation in the boron tetrahedron. The predicted 
area is 1560 (electrons) 2 which agrees to within 15 
percent with the observed area of 1750 (electrons) 2 . 

The 3.1 A peak appears to be a superposition of 
Ba-O, O-O, and Ba-B peaks. The Ba-0 distance 
predicted on the basis of ionic radii is 2.8 A. Con- 
sidering the distortion of the Ba cube, however, a 
larger distance is not unreasonable. The predicted 
area of the Ba-0 peak is 1698 (electrons) 2 . This 
figure is based on a value of 54 scattering electrons 
for barium. The 0-0 distance predicted is 3.1 A. 
This distance is the longer of the edges of the dis- 
torted barium cube The predicted area of this 0-0 



peak is 1040 (electrons) 2 . The predicted Ba-B dis- 
tance is one-half the Ba-Ba distance, or 3.4 A. The 
predicted area of this peak is only 173 (electrons) 2 . 
Thus, the peak makes no observable contribution to 
the radial distribution beyond its contribution to the 
area of the 3.1 A peak. The predicted total area of 
the 3.1 A peak is 2811 (electrons) 2 , in good agree- 
ment with the observed value of 2779 (electrons) 2 . 

The 4.3 A peak corresponds to the distance between 
a barium and the closest oxygens which are not near- 
est neighbors. These oxygens are the parts of the 
boron tetrahedron which are not involved in the edge 
sharing with the polyhedron of the barium under 
consideration . 

The 5.0 A peak corresponds to the distance across 
a face diagonal of the distorted barium cube. In 
addition, there should be a peak at 5.6 A correspond- 
ing to a body diagonal in the barium polyhedron. 
The area of this peak would be only 520 (electron ) 2 . 
It is not observed. 

Figure 7 shows the atomic radial distribution 
function obtained from the data. The most inter- 
esting feature of the peak is the clearly visible 6.8 A 
peak. Peaks at this distance 4 are not usually ob- 
served in glasses. It seems reasonable to assume 
either that the introduction of barium into the glass 
has increased the range of ordering or that because 
of barium's great scattering power it allowed the 
determination to be extended to this region. 
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Figure 7. Atomic radial distribution function, 47rR 2 p(R), 
calculated from the data of figure 8. 
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Because of the agreement between the position of 
peaks on the radial distribution and the distances 
predicted on the basis of structure type B, it seems 
likelv that the structure is type B and is definitely 
not type A. The good agreement between predicted 
and observed areas i\(\d* validity to this conclusion. 
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